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Abstract. Models describing transport and diffusion processes occurring 
along the edges of a graph and interlinked by its vertices have been 
recently receiving a considerable attention. In this paper we generalize 
such models and consider a network of transport or diffusion operators 
defined on one dimensional domains and connected through boundary 
conditions linking the end-points of these domains in an arbitrary way 
(not necessarily as the edges of a graph are connected). We prove the 
existence of Co-semigroups solving such problems and provide conditions 
fully characterizing when they are positive. 


1. Introduction 

Recently there has been an interest in dynamical problems on graphs, where 
some evolution operators, such as transport or diffusion, act on the edges 
of a graph and interact through its nodes. One can mention here quantum 
graphs, see e.g. [101 HH [Ml US] , diffusion on graphs in probabilistic context, 
unmaiM], transport problems, both linear and nonlinear, miiiiiTiiTiiMi 
HO] , migrations, [M], and several other applications discussed in e.g. [MllM]. 
In particular, the recent monograph |27] is a rich source of network models 
and methods. However, most of these works focus on a particular type of 
problems. For instance, in the quantum graph theory the main interest is to 
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determine whether the operators defined on the edges of a graph are self- 
adjoint and the work is confined to the Hilbert space setting. Most papers 
on the linear transport theory on graphs focus on long term dynamics of the 
flow. Papers such as [iniiiiiEg, motivated by probabilistic applications, look 
at Feller or Markov processes on graphs. 

The present paper, which provides the theoretical foundation for 0, is 
similar in spirit to [inille] in the sense that we prove the existence of strongly 
continuous semigroups in the space of continuous functions, as well as in the 
space of integrable functions, that solve the diffusion problem on a network. 
However, we extend the existing results of [iniiii] by considering processes 
that are more general than the diffusion along edges of a geometric graph 
with Robin boundary conditions at its vertices in the sense that we allow for 
a communication between domains that are not necessarily physically con¬ 
nected. In fact, the models we analyse can be also interpreted as diffusion 
on a hypergraph, , but we shall not pursue this line of research in this pa¬ 
per. For completeness, we also present similar results for transport problems, 
generalizing [5] in a similar way. 

To explain the idea of our extension, in the next section we consider two 
examples, see (gun]. 

1.1. Motivation 

First, let us introduce basic notation which will help to formulate the prob¬ 
lems and results. We will work in a finite dimensional space, say, R™. The 
boldface characters will typically denote vectors in R™, e.g. u = (ui, ..., Um)- 
We denote Ad = {1,..., m}. Further, for any Banach space W, we will use 
the notation X = X x ... x X, e.g. for X = Li(I),I = [0,1] we denote 

m times 

Li(/) = Li(J)x...xLi(/). 

'-V-" 

m times 

Let {A,D{)) be an operator in X. If A is a generator, we will denote 
by {e‘^}t>o the semigroup generated by A. 

1.1.1. Diffusion. We consider a finite graph without loops and isolated edges 
Q = (V, £) with, say, n vertices and m edges. On each edge there is a substance 
with density Uj,j G Ad, which diffuses along this edge according to 

dtUj — (7jdxx'^j 1 ^ f [) J ^ -'^7 (d) 

where Uj > 0 are constant diffusion coefficients, and can also enter the adja¬ 
cent edges. To simplify considerations, each edge is identified with the unit 
interval /. In the model of m, the particles can permeate between the edges 
across the vertices that join them according to a version of the Fick law. To 
write down its analytical form, first we note that, since diffusion does not 
have a preferred direction, we can assign the tail, or the left endpoint, (that 
is, 0) and the head, or the right endpoint (that is, 1) to the endpoints of the 
edge in an arbitrary way. Let li and be the rates at which the substance 
leaves through, respectively, the left and the right endpoints and Uk and 
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rik be the rates of at which it subsequently enters the edge Cfc. Then the Pick 
law at, respectively, the head and the tail of e^, gives 

-^xU^{l) = riUi{l) - y^JijUjjv), 

^ (^) 

where we have written Uj{v) as v may be either the tail or the head of the 
incident edge Cj. In particular, if there are no edges incident to the tail of e^, 
or there are no edges incident to the head of e^, then the Tick’s laws take the 
form 

dxUi{0) = kui{0), -da:Ut{l) = riu,{l), ( 3 ) 

respectively, where either coefficient on the right hand side can be 0. 

It is clear that if ^ 0, then lij = 0 and if lij ^ 0 then = 0 and 
thus we can define aw my- m matrix A by setting = 1 if either = 1 or 
lij = 1 and zero otherwise. Then A is the adjacency matrix of the line graph 
L{Q) of see e.g. [H]. However, it turns out that such a matrix is not easy 
to use and we see that introducing, for any i,j G Ai, 

kij = —lij if ?; = 0, kij = —lij if u = I, k^i = k, 
kij = rij if ■y = 0, kij = if u = I, = -U, (4) 

where v is either the tail or the head of the edge under consideration, the 
problem can be written as 


dtu{x,t) 

= Dda;xu{x,t), (a;,t) ejO, l[xR+, 


d^u{0,t) 

= K°°u(0,<)+K°iu(l,t), 

t > 0, 


dxu{l,t) 

= K^°u(0, <) + K^^u(l, t). 

t > 0, 


u(a;, 0) 

= u(x), xe]0,l[. 


(5) 


where u = (mi, ..., Um), D = diag{(Ti}i<i<m and u is the initial distribution. 

It turns out that there is no mathematical reason why the matrices K ‘^, 
wGH = {00,01,I0,ll},in © should be restricted to the matrices given by 
which, indeed, form a strictly smaller class, see [S]. In this paper we study 
the well-posedness of ([5]) for arbitrary matrices K‘^ in spaces C(/) and Li(/), 
extending and simplifying the results of [TOlIIS]. We also find necessary and 
sufficient conditions for the semigroup solving ([S]) to be positive. 

1.1.2. Transport problems. We consider a digraph G = (y{G), E{G)) = 
({ui,..., Vn}, {ei,..., Cm.}) with n vertices and m edges. We suppose that 
none of the vertices is isolated. As before, each edge is normalized so as 
to be identified with I with the head at 1 and the tail at 0. Following 
[3 m m m US], we consider a substance of density Uj{x,t) on the edge 
Bj, moving with speed Cj along this edge. The conservation of mass at each 
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vertex is expressed by the Kirchhoff law, 

m m 

i=i i=i 

where $“ = {4>ij)i<z<n,i<j<m and $+ = {4>tj)l<^<n,l<j<m are, respectively, 
the outgoing and incoming incidence matrices; that is, matrices with the entry 
(j)~j (resp. cj)fj) equal to 1 if there is edge ej outgoing from (res. incoming to) 
the vertex Ui, and zero otherwise. Note that due to definitions of the matrices 
and $+, the summation on the right hand side is over all incoming edges 
of the vertex Vi and on the left hand side over all outgoing edges of Vi . 

In [3] we considered a slightly more general model 

dtUj(x,t) + CjdxUj(x,t) =0, X G (0,1) , t >0, 

Uj{x, 0) = Uj{x), 

m 

t) = Wij (j)+^{jkCkUkil, t)), (7) 

fc=i 

where > 0 and > 0 are the absorption/amplification coefficients at, 
respectively, the head and the tail of ej. Here the matrix {wy }i<i<rt,i<j<m 
describes the distribution of the incoming flow at the vertex Vi into the edges 
outgoing from it; it is a column stochastic matrix, [53]. We denote C = 
diag{cj}i<j<m,5 = diag{^j}i<j<m and P = diag{ 7 j}i<j<m- 

It follows, [3], that if G has a sink, than there is no Co-semigroup solving 
0 . Hence we discard this case and then it can be proved, e.g. m Proposition 
3.1], that o can be written as an abstract Cauchy problem 

Ut = Au, u(0) = u, (8) 

in X = Li(/), where A is the realization of A = dia,g{—Cjdx}i<j<m on the 
domain 

D{A) = {u G Wliiy, u(0) = 5-iC-iirCu(l)}, (9) 

where B is the (transposed) adjacency matrix of the line graph of G. 

As with the diffusion problems, there is no mathematical reason to re¬ 
strict our analysis to the matrices of the form in the boundary 

conditions which, in fact, see |S|, form a strict subset of the set of all matrices. 
Thus, we consider the following generalization of (jHI, (Hj), 

= Au, u(0) = IKu(l), u(0) = u, (10) 

where K is an arbitrary matrix. 

Example 1.1. The main difference between m with arbitrary K and the 
model with K given in m is that in the former, the exchange of the sub¬ 
stance can occur instantaneously between any edges, while in the latter the 
edges must be physically connected by vertices for the exchange to take place. 
So, for instance, if there is a connection ei —>■ 62 and 62 —t 63 , then there is no 
connection ei —^ 63 . So, while in general m cannot model a flow in a physi¬ 
cal network, it can describe e.g. a mutation process. Indeed, let a population of 
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cells be divided into m subpopulations, with v(x) = (z;i(a;),... ,Vm(x)), where 
Vj{x),j S M.,x € [0,1], is the density of cells of age x whose genotype be¬ 
longs to a class j (for instance, having j copies of a gene of a particular type). 
We assume that cells of class j mature and divide upon reaching maturity at 
x = 1, with offspring, due to mutations, appearing in class i with probability 
kij, i G M. In such a case is a column stochastic matrix. We 

note that a particular case of this model is the discrete Rotenberg-Rubinov- 
Lebowitz model, [30]: where the cells are divided into classes according to their 
maturation velocity. 

A similar interpretation can be given to W where the variable x, instead 
of the age, denotes the size of the organism, e.g. [18j . 

Moreover, problems of the form m arise e.g. in queuing theory, [T3] . 


2. Well-posedness of the diffusion problem 

We shall consider solvability of (O in X = C(/) and X = Li(/). For technical 
reasons, we also shall need the solvability of (l5|) in Wi[(/). The norms in these 
spaces will be denoted, respectively, by || • ||oo, || • ||o: II • ||i- However, if it does 
not lead to any misunderstanding, we will use X to denote any of these spaces 
and II • II or || ■ ||x to denote the norm in X. Similarly, by ||| • |||x we denote 
the operator norm in the space of bounded linear operators from X to X. 

Our results are based on m and thus, introducing relevant spaces and 
operators, we try to keep notation consistent with op.cit. First, consider 


X 3 u ^ Lu = (7 o5,,u,7i5,,u) g Y = R"* x R"*, (11) 


where 7^, z = 0,1, is the trace operator at x = i (taking the value at a: = z if 
X = C(/)). The domains of L are D{L) = C^I) if X = C(/) and D{L) = 
W^(/) in two other cases. Then we define the operator 

X 3 u ^ chu = K(7oU,7iu) = ( j|io In ) ( ^ x R-. 


Further, let us denote 

<I>*u = K*(7 oU, 7 iu) 


j^OOT _J^10T W 70U \ 

y V 71" / ’ 


( 12 ) 


where denotes the transpose of K. Clearly ($*)* = $. 

Let A denote the differential expression Au := I[!)9a;a:U. Then we define 
the operators A^, a = 00, 0,1 by the restriction of A to the domains 


D(A^) = {uG C2(/); Lu = $u}, 

D{Al) = {uG W2(/); Lu = $u}, 

D{Al) = {uG W3(/); Lu = $u}, (13) 


respectively. As before, we drop the indices from the notation if it will not 
lead to any misunderstanding. 
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2.1. Basic estimates in the scalar case 

Consider the general resolvent equation for ([S]) 

Au - = f, a:G]0,1[. (14) 


Since without the boundary conditions the system is uncoupled, its solution 
u is given by 

u(x) = E(—/ra;)Ci + E(/rx)C2 + U/x(a;) (15) 

where E(±/ra:) = diag{e=*=''*“}i<i<m, 0 ^ Xjoi = fif = |A/tTi|e*® with Kyui > 
0, U^(x) = (U^i (x),..., (x)) with 

1 

UM = 7^ [ e-'^^l“-^l/.(s)ds, (16) 

0 

and the vector constants Ci and C2 are determined by the boundary condi¬ 
tions. Further, denote 

= {A = |A|e*® G C; |A| > a, \e\ < 9o < tt} 


and . 

The starting point are well-known estimates for the scalar Dirichlet and 
Neumann problems, e.g. m- We briefly recall them here in a slightly more 
precise form, similarly to [S]. In the scalar case we can use cr = 1 as will 
become clear when the estimate is derived. It follows that for the Dirichlet 
problem we have 


1/^(1)-t/^(0)e^ 

1 g/i _ g-M 

and for the Neumann problem 

^ U^{0)e^ + C/41) 


Further, m, 
\\UJoo < 


gM - g-M 


< 


^ _ U^{0)e-^ - U^{1) 


C^ = 


g/i _ g-/i 


U^jl) + l/^(0)e-^ 

2 g/i _ g-^ 


|/r|5R/r |A|cos0/2 

and, for z = 0,1, 


IlCfJlo < 


< 


|A|cos0/2 




(17) 


(18) 


2|/i|3?/r ’ ^ 2|^| ■ 

The following result plays an essential role in deriving precise estimates. 


(19) 


g/i _ g-/i 




E' 

n—O 


,-2n/i 








n—0 


I _ ' 


( 20 ) 
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Then, for either Dirichlet or Neumann problem in C{I), for A S Egp for any 
9o < TT, we have for ui = D,N 




M b 


< 


< 

ll/lloo ( 

A cos do/2 \ 

< 

2||/l|o 

A cos do/2 

in Li 

(I), we have 

< 

1 _ 

< 

Il/llo ( 

A cos do/2 \ 

< 

2||/l|o 


|A| cos do/2 
1 


- 1 


. — e 


-2SR/i 


+ 1 


( 22 ) 




+ 1 


I /\ I COS t/Q / z 

Consider now Aq, see (USD; that is, the operator corresponding to the Neu¬ 
mann boundary conditions in one dimension. We have 


Proposition 2.1. Aq generates an analytic semigroup in W/(/) with the re¬ 
solvent satisfying the estimate 

||i?(A, Al)f\\, < A e (23) 

IAI cos (yQ/ 2 

Proof. Consider the resolvent equation 

Am - da;xU = /, 9a:M(0) = dxu{l) = 0. 

If / S Wf{I), then M G Wi{I) and we can differentiate the differential 
equation getting, for v := dxU, 

Am - dxxV = dxf, m(0) = m(1) = 0; (24) 


that is, V satisfies the resolvent equation for the Dirichlet problem. Hence 

2||/||o , 2||a,/||o 2||/||i 


|l = ||u||o -I- ||5a;M||o < 


-b 


|A| cos do/2 |A| cos do/2 |A| cos do/2’ 


□ 


Since Am — adxxU = / is equivalent to cr ^Am — dxxU = a ^/, we see 
that the estimates above are independent of cr. 

2.2. Solvability of ([5]| 

The ideas in this section are based on [inmi but the analysis is simplified by 
using the analyticity of the semigroup and the application of [m Theorem 
2.4]. 
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Since in the vector case and Neumann boundary conditions the resolvent 
equation decouples, we obtain 

||i?(A, Ao)/||x < A e (25) 

IAI cos OqI2 

for any 9o < tt and X = C(/),Li(/), W)(J). Hence, in particular Aq gener¬ 
ates an analytic semigroup in X. 

We begin with a straightforward consequence of Ref. m- 

Theorem 2.2. Let X be either C(/), or Wi[(/). Then the operator A$ gen¬ 
erates an analytic semigroup in X with the resolvent satisfying the estimate 

|||i?(A,A,s)|||x <2|||R(A,Ao)|||x, A G (26) 

for some a > 0. 

Proof. The result follows directly from [m Theorem 2.4], since Aq generates 
an analytic semigroup in X, L is an unbounded operator on X which is, 
however, bounded as an operator from D{A) to Y, and it is a surjection (the 
right inverse in each case is given by [^“^(yo, yi)](a:) = ^(yi - yo)a;^ - yoa;. 
Furthermore, $ is a bounded operator on X (if X = W](/) this follows since 
W](/) functions are absolutely continuous on I). Hence, the assumptions 
[T71 (1.13)] are satisfied and the theorem follows from [TT] Theorem 2.4]. □ 


In X = Li(/) the situation is more complicated as $ is not bounded on 
X. First we show that R{X, A^) extends to a resolvent on Li(/). 


Lemma 2.3. We have 



= A 


0 

$ 


and the resolvent set of satisfies p(A^) n IR.+ ^ 0. 


Proof. Let D{A\,) 9 u„ —>• u G Li(/) and A^u„ = Au„ —>• v G Li(/). This 
shows that u„ —>■ it in Wf (/) and thus v = Au. Since taking the trace of a 
Wf(/) function and of its derivative is continuous in Wf(/), the boundary 
values of u„ and 9xU„ are preserved in the limit and thus u G D{A^). Hence 



C At 


On the other hand, let u G Il(At■ Then v = u — f GkFi(^)) where 

f(a:) = (u(0)-I-u'(0) - 2u(l)-b u'(l))a;^ 

— (2u(0) -b 2u'(0) — 3u(l) -b u'(l))a;^ -b u'(0)x -b u(0) G Il(At ■ 


Thus there is a sequence (h„)„gN C C^{I) converging to v in Wf(/) and 
hence u is the limit in WKl) of functions u„ = h„ -b f G D{A]^). Since the 
convergence in implies the convergence of both u„ and Au„ in Li(/) 

we see that also 


Ai 


■Li(H 


D A! 


$• 
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Finally, we see that, as in the scalar case HMD, the solvability of the resolvent 
problem for m is equivalent to solvability of the linear system 

-MC 1 +MC 2 = ]K°°(Ci + C 2 )+K°i(e"''Ci + e'"C 2 ) 
-MU^(O) + K°°U^(0) + K°iU;,(l), 

-M]E(-/x)Ci+ME(/ x)C 2 = Ki°(Ci + C 2 )+ K^^(E(-/x)Ci+E(/ x)C2) 

+MU^(l)+]K°iU^(0)+KiiU;,(l), (27) 

for the vectors Ci, C 2 , where M = diag{/x}. Once the constants are found, the 
solution is given by the vector version of m and thus belongs to W^(/) C 
Li(/). The system above is solvable provided its determinant is different from 
zero. However, the determinant is clearly an entire function in ^ and thus 
can have only isolated zeros in C. Hence, there must be positive values of 
AGp(AO). □ 

Theorem 2.4. The operator generates an analytic semigroup in Li(/). 

Proof. We use the formula from Lemma 1.4 of Ref. stating that 

R(A,Ai) = (/-LA$)-iR(A,A]v) (28) 

in X = W)(/) as (/ — La 4*)“^ is an isomorphism of X for large |A|, where 
Lx = (L|/fer(A-A))~^, A G p{A]^), See the proof of [171 Theorem 2.4]. 

However, i?(A, A]^) extends by density to R{X,A^) (the resolvent on 
Li(J)) which is a bounded linear operator from Li(/) to I?(A]y) C W 2 (/) 
which is continuously embedded in W](J) and thus i?(A, A^) extends to a 
bounded linear operator, say i?$(A), on Li(/). Since R{X, A^) is a resolvent 
on W](/), we obtain, by density, that i?$(A) is a pseudoresolvent on Li(/). 
Furthermore, C“(]0,1[) C Il(A^), thus it is also a subset of the range of 
7?$ (A) and therefore the range is dense in Li(/). Also, since the range of 
R{X,A%) is in Wf(/) C W](/), there is no need to extend (/ — Lx^)~^. 
Hence i?$(A) is a one-to-one operator as a composition of two injective op¬ 
erators. Thus, by Proposition HI.4.6 in Ref. [15], 7?$ (A) is the resolvent of a 
densely defined operator, say A$. Clearly, A$ is a closed extension of A^ and 

thus A$ D A^ = A^ ^ . From Lemma [2j3l there is A G p(A$) fl p(A^), 

thus 7?(A, A$) D 7?(A, A^) and hence A$ = A^ and, consequently, 7?(A, A^) 
is defined in some sector. To prove that is a sectorial operator, we use the 
idea of m but in a somewhat simpler way. We consider the operator A^,. 
Let Rf denote the adjoint to 7?(A,Agl) for A G p(A“.), which acts in the 
space of signed (vector) Borel measures on J, [9]. R{X, A“.) is given by (fTKll 
with Ci,C 2 given by (1271) with the matrices K“ replaced by corresponding 
matrices in Hence, apart from U^(x), 77(A, A“.) is a composition of 
an algebraic operator coming from inverting the matrix in (HZl) with a vector 
of functionals acting on f. Thus, if f G Li(7) is the density of an absolutely 
continuous measure, a standard calculation shows that 

7?ff = 7?(A,A°)f 
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From the definition of the norm of a signed Borel measure, if the latter is 
absolutely continuous, its norm is equal to the Li norm of its density. Since 
taking the adjoint preserves the norm of the operator, we obtain 

|||i?(A,A°)|||L,(,) = |||i?(A,A-)|||c(/). (29) 

Hence, by Theorem [221 is sectorial and, being densely defined, it gener¬ 
ates an analytic semigroup on Li(J). □ 

2.3. Positivity of the semigroup 

Let us recall that for an element u of a Banach lattice, we write u > 0 if 
0 ^ u > 0. We have the following result 

Theorem 2.5. The semigroup {e*'''"* }t>o is resolvent positive if and only if 
— K™, are nonnegative off— diagonal 

—IK°^, are nonnegative. (30) 


Proof. To prove the result we use Theorem B-II.1.6 in [5H], which states that 
if an operator A on C{K) (where K is compact) generates a semigroup, then 
the semigroup is positive (or, equivalently, A is resolvent positive) if and only 
if A satisfies positive minimum principle: for every 0 < / < D{A) and x G A, 
if f{x) = 0, then {Af){x) > 0. However, to be rephrase the problem at hand 
in the language of the positive maximum principle, C{K) must to be the space 
of real valued functions. To achieve this, here we identify C(/) with the space 
(7(1), where I = [Oi, li]U.. .U[0m, Im]; that is, instead of considering a vector 
function on I we consider a scalar function on a disconnected compact space 
composed of m disjoint closed intervals. In particular, each edge e^-, j G Ad is 
identified with the closed interval [Oj, 1^]. Then A“ will be changed to A’^ 
which is the restriction of 

{Au){x) = X[o,,i^]{x)crjd^^u{x) 

j&M 

will be the space of twice differentiable functions on Inti =]0i, li[U ... U]0m, Im 
differentiable on I and satisfying 

m m 

dxu{Qj) = ^/c°)!u(0fc)-t-^fc°^u(lfc), 

k^l k^l 

m m 

dxu{lj) = ^fcj^u(Ofc) +^fcj^M(lfe). 

k—1 k—1 

Thus, if 0 < u G D{A) takes value 0 at some a; G I, then either x G Inti and 
then classically dxxu{x) > 0, or a: = 0^ or a: = Ij for some j G Ad. If a: = Oj, 
then dxu(Oj) > 0. If dxu(Oj) = 0, then dxxu{Qj) > 0 follows as in Example 
B-H.I.24 of Ref. [28] (or simply by noting that the even extension to [—1, 0] 
gives a function on [—1,1] with minimum at a: = 0). If dxu{{)j) < 0 then, 
since u(0j) = 0, such a function cannot be nonnegative on [Oj, Ij]. Analogous 
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considerations hold if m(1j) = 0. Therefore the positive minimum principle is 
satisfied and yields the positivity of }t>o and thus of }t>o. 

To prove the converse we introduce the following notation. For any 
a := (a°, a^) > 0 with = 0 for some i G A4,j = 0,1, we have 


Accordingly, we denote 


-K°° -KOI \ 

K^o K“ J • 


(31) 


(Sa): = (-1)’'+! ^ r = 0,l,seAd. 

i=o.i \jeM J 

Let us assume that (15(11) is not satisfied. Then there is a non-diagonal element 
of S which is strictly negative. Suppose (—1)’'+^A:[J < 0 for some i j and 
r, s = 0,1, and consider a vector a with 


= 1, aC = 0, = 5 > 0, f = 0,1, I ^ j it t = s and I ^ i it t = s. 

(32) 

Then for r = s we obtain 


(Sa)C = (-l)’'+M^-bd ^ 

where f = 0 if r = 1 and f = 1 if r = 0, and for r ^ s 



< 0 


(33) 


(S«)r = (-l)’'+M^-(-J ^ k:[+ fc? <0 (34) 

\ / J 

for sufficiently small d > 0. We shall prove that there exists a function 0 < it € 
C^)!) satisfying u{ri) = a\ and {—lY~^^dxu{ri) = (Sq:)C which additionally 
satisfies dxxu{ri) < 0. 

For a given constants aY Pi = (Sec))’,!’ = 0,1, we consider auxiliary 
functions /[{x) = Plx{x — l)-\-a\. We have /[{r) = al,dxfi{r) = (—1)’'+^/3[ 
and dxxfi = 2/3[. We observe that as long as > 0, there is a one-sided 
interval ((0,Wi) if r = 0 and {oJi,l) it r = 1), where /f > 0 irrespective of 
the sign of /3[. On the other hand, if = 0, for local nonnegativity we need 
PI < 0. Now let (/) be a nonnegative (7°° function which is 1 on [—a, a] and 0 
outside [—2a, 2a] where 0 < 2a < min^g^v^ Wi and define 

u{x) = (l){x)f°{x) + Pit - x)fj{x), X G [Oj, lj],j G M. 

is a C'°°(I) function that satisfies: 

1. It > 0; 

2. m(0j) = a° and dxu{0j) = dxfjiOj) = -P° = 

3. u(lj) = aj and dxu(lj) = dxfj(lj) = P] = (5a)], 

so that 0 < It G D{A). Recalling that we assumed (—l)’'+^fc’J < 0 for some 
i Y j and r,s = 0,1, we consider coefficients {a]}t=o,i,je 7 W satisfying (1321) . 
Then we have u{ri) = 0. On the other hand, dxxu{ri) = 2/3® = 2(Sa)C < 0 
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by (1551) or (1551) . Thus, there is a nonnegative element u £ D{A) for which 
dxxU < 0 at a point where the global minimum of zero is attained. □ 

We can use this result to prove an analogous result in Li(/). 

Corollary 2.6. The operator generates a positive semigroup if and only if 
the assumptions of Theorem \2 .51 are satisfied. 

Proof. In one direction the result immediately follows by density of C(/) in 
L^(J). Conversely, if {e‘"^'<'}t>o > 0 then, in particular, for any 0 < u £C(/) 
we have e*^'^u = e*^^ u > 0. □ 

3. Solvability of ([8]) 

We return to the problem Q where, we emphasize, K is an arbitrary matrix. 
In this section we restrict our attention to X = Li(/) as in C(/) the operator 
A is not densely defined. Let us recall that A is the realization of A = 
dda,g{—Cjdx}i<j<m on the domain D[A) = {u £ W}(/); u(0) = Ku(l)}. 

The following theorem for 0 has been proved in [5] (see also [4]) but the 
proof for any nonnegative matrix K is practically the same. Here we extend 
this proof to an arbitrary K. However, for the proof in the general case we 
need to provide basic steps of the proof for nonnegative matrices. 

Theorem 3.1. The operator (A, Z1(A)) generates a Co-semigroup on Li(/). 
The semigroup is positive of and only z/K > 0. 

Proof. Clearly, Cg°(]0,1[) C D{A) and hence D{A) is dense in X. To find 
the resolvent of A, the first step is to solve 

Xuj + CjdxUj = fj, j = I,..., m, X £]0,1[, (35) 

with u £ D{A). Integrating, we find 

X 

u(a:) = EA(a;)v + C~^ J Ea(x — s)i{s)ds, (36) 

0 

where v = (rii ,...,?;„) is an arbitrary vector and E^ (s) = diag |e £ ^ 

I J l^j^m 

To determine v so that u £ D{A), we use the boundary conditions. At x = 1 
and at X = 0 we obtain, respectively 

1 

u(l) = Ea(1)v + J Ea(1 — s)i{s)ds, u(0) = v. 

0 

Then, using the boundary condition u(0) = Ku(l) we obtain 

1 

(I - KEa(1))v = KC-i J Ea(1 - s)i{s)ds. (37) 

0 
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Since the norm of Ea( 1) can be made as small as one wishes by taking large 
A, we see that v is uniquely defined by the Neumann series provided A is 
sufficiently large and hence the resolvent of A exists. 

Let us first consider K > 0. Then the Neumann series expansion ensures 
that A is a resolvent positive operator and hence we can consider only f > 0. 
Adding together the rows in (1371) we obtain 


i=i 


= H 

i=i 


Vj + 


E- 

i=i 




(38) 


771 


where Hj = ^ kij. Renorming X with the norm ||u||c 
2 — 1 

and using dSH) and we obtain 


E m 

.7 = 1 


^-1| 




= iE 




i=i 


i=i 




1 




We consider three cases. 

(a) Kj < 1 for j G M. Then ||Ea(—1)K|| < 1 and thus v, and hence R{X,A), 
are defined and positive for any A > 0. Further, dropping the first two terms 
in (IMl) we get 


u 


i=i 


If 1 

- J fj{s)ds = j\\f\\c, A>0. 


Hence (A, D{A)) generates a positive semigroup of contractions in (X, || • ||c). 

(b) Kj > 1 for j G Ad. Then (IM|) implies that for some A > 0 and c = 1/A we 
have 

||i?(A,A)f|U>c||f|U 

and, by density of D{A), the application of the Arendt-Batty-Robinson the¬ 
orem mm, gives the existence of a positive semigroup generated by A in 
(X, II • lie). Since || • ||c is equivalent to || • ||x, A generates a positive semigroup 

in X. 

(c) Kj < 1 for j G h and Kj > 1 for j G h, where Ii Ci I 2 =0 and 
/i U /2 = {1,. •. ,m}. Let L = {lij)i<ij<m, where kj = kij for j G h and 
kj = 1 for j G Ii. Denoting by Al the operator given by the differential 
expression A restricted to D{Al) = {u G W((/); u(0) = Lu(l)} we see, by 
(|3^ . that 

0 < R(A, A) < i?(A,AL) (39) 

for any A G /^(Al). By item (b), Al generates a positive Co-semigroup and 
thus satisfies the Hille-Yosida estimates. Since clearly (15^ yields i?^(A, A) < 
i?^(A, A]l) for any fc G N, for some w > 0 and M > 1 we have 

||i?'=(A,A)|| < ||R'=(A,Al)|| <M(A-w)-^ A>a;,fcGN, 


and hence we obtain the generation of a semigroup by A. 
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Assume now that K is arbitrary. The analysis up to (1371) remains valid. 
Then (IS71) can be expanded as 


OO p 

= ^(KEa(1))”KC-i / Ea( 1 - s)f(s)ds 

n=Q i 


(40) 


and hence 


OO " ^ 

u(a;) =EA(a:)^(]KEA(l))”]KC-^ / EA(l-s)f(s)ds+C-^ / EA(a:-s)f(s)ds. 

n=0 { { 

(41) 

Denoting now |u| = (|mi|, ..., |m„|) and |K| = {\kij\)i<ij<rn, and using the 
fact that only K may have non positive entries, we find 

OO \ 

|u(a;)| < Ea(x)^(|K|Ea( 1)|)”|K|C-Y EA(l-s)|f|(s)ds 

n=0 { 

X 

+ iC-^j¥.x{x-s)\i\{s)ds = i?(A,A 
0 

where A|k| denotes the expression A restricted to D(A|]k|) = {u G W}(/); u(0) 
|A|u(l)}. So, we can write 


|i?(A,AK)f| <i7(A,A|K|)|f| (42) 

and, iterating, 

|i?(A,AK)'=f| <i?(A,A|K|)'=|f|. 

Using the fact that taking the modulus does not change the norm, we find 

P(A,AK)'=f||<P(A,A|K|)'=|f|||<y^||f||. 

with M and uj following from the Hille-Yosida estimates for A|]k|. 

The fact that K > 0 yields the positivity of the semigroup follows from 
the first part of the proof. To prove the converse, let kij < 0 for some i,j and 
consider the initial condition f(x) = (/i(x),..., /m(x)) with fk = 0 for k ^ j 
and fj € (^^^([0,1]) with fj(0) = fj(l) = 0, so that f G D{A), and fj{x) > 0 
for 0 < X < 1. Then, at least for t < mini<j<m{l/cj}, Ui satisfies 

— CiOxUi , Ui (x, 0) — 0, Ui (0, t) — kij fj (1 Cjtf 


that is, 


(x,t) = kijfj 1 + - Cjt 


Ci 


X 

t > — 

Ci 


and we see that the solution is negative for such t. This ends the proof. □ 


We conclude the paper by noting an interesting corollary which also uses 
monotonicity properties of the problem and which is important in asymptotic 
analysis, see [B]. For given velocities C = (ci,... ,Cm), let c = minjg^{cj} 
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and Cmin = {c,..., c}. For a particular velocity matrix C, let Ak,c denote 

m times _ 

the generator of the semigroup solving (US. 

Corollary 3.1. There holds 

||gtAK,C|| < ||gtA|K|,C„i„||^ 

Proof. We have by (15^ and (liOl) that 

||i?(A,AK.c)f|| < ||i?(A,A|K|,c„,J|f|||, 

hence the result follows from [29l Theorem 1.8.3]. □ 
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